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Abstract 



We present a study of Goldstone boson scattering based on the flavor SU(3) chi- 
ral Lagrangian formulated with vector mesons in the tensor field representation. A 
coupled-channel computation is confronted with the empirical s- and p-wave phase 
shifts, where good agreement with the data set is obtained up to about 1.2 GeV. 
There are two relevant free parameters only, the chiral limit value of the pion decay 
constant and the coupling constant characterizing the decay of the rho meson into a 
pair of pions. We apply a recently suggested approach that implements constraints 
from micro- causality and coupled- channel unitarity. Generalized potentials are ob- 
C<*) tained from the chiral Lagrangian and are expanded in terms of suitably constructed 

conformal variables. The partial-wave scattering amplitudes are defined as solutions 
of non-linear integral equations that are solved by means of an N/D ansatz. 

o 

>■ 1 Introduction 
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The study of Goldstone boson interactions is a time honored challenge in 
hadron physics [1,2 3,4 5,6,7,8]. In the last decade there was significant progress 
in a profound description and understanding of the low-energy regime based 
on the chiral Lagrangian. In the close-to-threshold region chiral perturbation 
theory (%PT) is applicable and leads to controlled results that are consistent 
with the empirical scattering data PllU|lll|ll2|[T3"] . Extensions of xPT, that im- 
plement coupled- channel unitarity by means of partial summation techniques, 
aim at a description of the scattering data at somewhat larger energies where 
meson resonances play an important role [l4"fl5Jll6JfT7"] . Most remarkable is the 
systematic computation at the one-loop level where the inverse of the partial- 
wave amplitudes is expanded in chiral powers [17J. Accurate results for the s- 
and p-wave phase shifts up to about y/s ~ 1.2 GeV were obtained by adjusting 
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the Q 4 counter terms of the chiral Lagrangian. This implies that the lowest 
scalar and vector mesons with J p = + and 1~ are properly described by such 
an approach. It has been argued by detailed large- N c scaling studies that the 
nature of the light scalar and vector mesons is distinct [18]. While the scalar 
mesons may be generated by the leading order terms of the chiral Lagrangian, 
the vector mesons are a consequence of the subleading counter terms. 

The purpose of this Letter is to further explore the dynamic role of light 
vector mesons in a chiral Lagrangian. The motivation of considering the light 
vector mesons as explicit degrees of freedom is twofold. First, we recall that 
it was shown previously that the size of the symmetry preserving Q A counter 
terms may be estimated by a resonance saturation mechanism where the light 
vector mesons appear to dominate [TH]. Second, the light vector mesons play 
a particular role in the hadrogenesis conjecture [1812 0]. Together with the 
Goldstone bosons, they are identified to be the relevant degrees of freedom that 
are expected to generate the meson spectrum. For instance it was shown that 
the leading chiral interaction of Goldstone bosons with the light vector mesons 
generates an axial-vector meson spectrum that is quite close to the empirical 
one [20] • From this point of view the axial- vector mesons are similar in nature 
as the scalar mesons. We consider a chiral Lagrangian with explicit vector 
meson fields as a promising starting point to compute the meson spectrum. 

While it is straightforward to write down a chiral Lagrangian with explicit 
vector mesons it is an open issue how to order the various terms according 
to some counting scheme. In a recent work a counting scheme was suggested 
based on large- iV c arguments and the dynamical assumption of hadrogenesis 
[18]. So far it has been tested successfully against two and three body decays 
systematically at the tree level [T8"]|2~T] . In this work we scrutinize its implica- 
tion in the scattering of Goldstone bosons. At leading order the vector meson 
exchange processes should be considered already, rather than at a subleading 
order, as is implicit in the standard xPT approach. Based on this assumption 
we perform a coupled-channel computation, where we apply a novel unitariza- 
tion scheme developed by one of the authors recently [221123] . 




Though the inverse amplitude method applied with great success in previ- 
ous works obtains an inverse amplitude that can be represented in terms of 
dispersion-integral representation there is a subtle limitation of this approach, 
if applied to a coupled-channel situation. In the inverse amplitude method or 
any algebraic approach the partial-wave scattering amplitudes have unphysi- 
cal left-hand branch points. This holds at any finite truncation. The source of 
these is easily understood. The locations of the left-hand branch points of a 
partial-wave scattering amplitude depend decisively on the channel. There is 
no universal branch point that characterizes all coupled-channel amplitudes. 
In any algebraic method there is a determinant of functions with left-hand 
branch points at different locations. Via the determinant a branch point of a 
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given channel is transported into any other channel as long as the transition 
potential is non- vanishing. In general this leads to the presence of all possible 
left-hand branch cuts in all channels. Though this is unphysical it does not nec- 
essarily always lead to numerically significant effects in the physical region. 
If all considered branch points are below the smallest considered threshold 
typically the presence of unphysical branch points are not really problematic. 
However, once a left-hand branch point of a heavy channel is located right 
to the threshold of a lighter channel the limitations of an algebraic approach 
start to be visible. 

The main objective of the novel scheme [221123] to be applied in this work is 
a controlled realization of the causality and unitarity condition in a perturba- 
tive application of the chiral Lagrangian. The starting points are partial-wave 
dispersion relations. A generalized potential is constructed from the chiral La- 
grangian in the subthreshold region and analytically extrapolated to higher 
energies. The partial-wave scattering amplitudes are obtained as solutions of 
non-linear integral equations. 

In our scheme the empirical s- and p-wave phase shifts are recovered up to 
about 1.2 GeV. There are two relevant free parameters only, the chiral limit 
value of the pion decay constant and the coupling constant characterizing the 
decay of the rho meson into a pair of pions. 



2 Chiral Lagrangian with vector meson fields 



We recall the leading order chiral Lagrangian for the Goldstone bosons and 
vector mesons. The vector mesons are represented in terms of anti-symmetric 
fields V^y = —V UfM . It has been demonstrated first in [19], that vector mesons 
interpolated in terms of anti-symmetric fields provide a transparent satura- 
tion model for the Q counter terms of the chiral Lagrangian in its Goldstone 
boson sector. The chiral Lagrangian is composed out of the building blocks, 
Ufj,, V^v, x±, that transform identically under chiral rotations (see e.g. [24J). 
We do not consider electromagnetic interactions here. Using covariant deriva- 
tives, _D M , the transparent transformation properties are maintained for terms 
involving derivatives, 
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where the pseudo-scalar meson octet field $ is normalized as in [18]. There 
are two types of terms which one has to consider. The ones that break the 
chiral SU(3) symmetry explicitly involve the field combinations x± an d are 
proportional to the quark-mass matrix of QCD. At leading order the diagonal 
matrix xo can De expressed in terms of the pion and kaon masses as indicated 
in ([2]). Isospin breaking effects are neglected. We consider further constraints 
from QCD as they arise in the limit of a large number of colors N c |18j . This 
implies that interaction terms that involve a double trace in flavor space can 
be neglected. The latter are suppressed by 1/N C as compared to single-flavor 
trace interactions. To order Q 2 the following terms were proposed in [TS] to 
be relevant in our case 



£ = / 2 tr{f^f/t + _ Ur\{D»V, a ){D u Vn} (3) 

+ g m\- tr {y^ V^} + h D tr |^ V, v X+ ) + i tr \u, U u } 

where / ~ f n may be identified with the pion-decay constant, f n = 92.4 
MeV, at leading order. A precise determination of / requires a chiral SU(3) 
extrapolation of some data set. In |25j the value / ~ 90 MeV obtained from 
a detailed study of pion- and kaon-nucleon scattering data. This value was 
used consistently in various applications of chiral Lagrangians to meson and 
baryon resonance physics |20|I27|28] . We discriminate the mass of the light 
vector mesons m x - and the scale parameter my. From [T8"f2"T] we recall the 
value for the parameter, hp, as determined by two and three-body decay 
properties of the light vector mesons 

hp = 0.29 ± 0.03 for m v = 776 MeV (4) 

at tree-level computations. A first coupled-channel computation based on (j3]), 
was presented in [20] where the scattering of Goldstone bosons off the vec- 
tor mesons was considered. The axial-vector meson spectrum was generated 
dynamically. 



3 Coupled-channel dynamics 



Given the Lagrangian ([3]) we compute the two-body coupled-channel scatter- 
ing amplitudes. At tree-level the scattering amplitude takes the generic form 



4 



T(s,t,u) 



a 



Ct C u 2 Ck 2 

s + TTT^ t + — — - u + — — - 77* + ~ ~ 777 ; 



I2f 2 Yip 12p 



12/ 



2 i" 



12/ 



2 



(5) 



(m y /ip) 2 C ; (a;) 



s—ch 



x=0 128P s-m 

8 

E 



x=0 



128/ 4 * — 171= 



x=0 



128/ 4 M-mJ 



(s + m 2 — m 2 ) (s + m 2 . — m\) + 2 ( t — m\ — m 2 ,) s 
( t + m 2 — m 2 )(t + m x — m 1 ) + 2 (s — m 2 — mj t 
(u + — m 2 )(-u + — th 2 ) + 2 (£ — m 2 — m 2 ) w 



where the coefficients C... depend on the channel, in particular on isospin (/) 
and strangeness (S). We refrain from detailing the various coefficients since 
they can easily be calculated and mostly have been given previously [TT] . 
The sums in ^ extend over the members of the vector meson nonet, where 
m x denotes their respective masses. The sum of the Mandelstam variables 
s + t + u — Tn\ + ?77,2 m i m 2 i s the sum of the squared meson masses, where 
we use a bar for the final state masses. Partial- wave amplitudes are introduced 
by an average 



T J (s) = f +1 ^ (^—) " T (s, t, u) Pj(cos 6) , (6) 

J — 1 \ Pcm Pcm ) 

over the center-of-mass scattering angle 9. The relative momenta of the initial 
and final state are denoted by p cm and p cm . The total angular momentum is 
J. 

The coupled-channel partial- wave amplitudes T^ b (s) are determined as solu- 
tions of the non-linear integral equation 



rpj / \ _ ttJ ( a \ I f°° ^ S S L^M Tg C (s) Pc d (s) Tfa (s) 



with the phase-space matrix 



The matching scale n 2 M is identified with the smallest two-body threshold 
accessible in a sector with isospin and strangeness (I, S). 

Given a generalized potential U^ b (s) the non-linear integral equation implies 
a partial- wave scattering amplitude T^ b (s) that complies with the coupled- 
channel unitarity request and the microcausality condition. On the other hand 
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the crossing symmetry constraint is not automatically satisfied [29]. If crossing 
symmetry would be implemented exactly, a crossed reaction can be computed 
from the direct reaction amplitude. Necessarily, the crossed reaction would be 
determined by the direct reaction amplitude evaluated at energies below the 
s-channel threshold. Now, imagine that in a first step we have some means to 
approximate the generalized potential systematically and very accurately at 
energies above the s-channel threshold only. We may use this potential and 
determine via ^ the scattering amplitude in the physical region. In this case 
we are not able to compute the crossed reaction, since our approximation for 
the generalized potential is valid above threshold only, by assumption. The 
crossed reaction may be computed by setting up the analogous non-linear 
integral equation ^ for the crossed reaction. It involves potentials distinct 
from the potentials of the direct reaction. But where is the crossing symmetry 
constraint? Our argument illustrates that crossing symmetry is a constraint 
that affects dominantly amplitudes at subthreshold energies. In a second step 
we now imagine that we have some means to approximate the generalized 
potential systematically at energies \fs > A , i.e. we include a small sub- 
threshold region. In this case crossing symmetry does provide a constraint. 
The crossed reaction amplitude can be computed from the direct amplitudes, 
however, only in a specific subthreshold region. The coincidence of the cross- 
ing transformed amplitudes and the amplitudes from the crossed reaction in 
that specific subthreshold region defines the desired constraint. Note that a 
non-empty coincidence region requires A to be sufficiently distinct from the 
s-channel unitarity branch point. In our scheme a measure for the amount of 
crossing-symmetry violation is the strength of non-perturbative contributions 
to the subthreshold scattering amplitudes. 

The non-linear integral equation ([7]) does not necessarily have a solution or 
it can have several. A unique solution is singled out by the condition that 
in a small vicinity of the matching point the results of chiral pertur- 
bation theory are recovered. This implies a significant suppression of possible 
crossing-symmetry violating contributions. From the form of ^ it follows that 
the existence of a solution requires the generalized potential to be bounded 
asymptotically, modulo some possibly logarithm terms. Thus, an evaluation 
of U^ b (s) in xPT is futile. Any finite order truncation leads to unbounded 
potentials, characterized by an asymptotic growth in some power of s. How- 
ever, as was pointed out in [22,23], the generalized potential U^ b (s) can be 
reconstructed unambiguously in terms of its derivatives at a point y? ah E that 
lies within its analyticity domain and where the results of %PT are reliable. A 
solution of ([7]) requires the knowledge of the generalized potential, U^ b (s), for 
real energies larger than the maximum of the initial and final thresholds only. 
Following [22] we identify with the mean of initial and final thresholds. A 
Taylor expansion of U(s) around /i E has a small convergence radius, that is 
determined by the distance to the closest left-hand cut. In our case the clos- 
est left-hand cuts reflect the two-body t- and u-channel unitarity condition. 
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In order to extend the convergence up to some cutoff scale A s , we apply the 
conformal map, that was constructed in [22]- Explicitly, it is given by 



aS) "(a-2&)(A2- s ) 2 + f a ~ (Af^lF' " (A] - Ag)^ ' (9) 

where the parameter A is identified unambiguously such that the mapping do- 
main of the conformal map touches the closest left-hand branch point. Within 
this domain the generalized potential can be reconstructed in terms of its 
derivative at the expansion point fi 2 E . It holds 



tf00 = £c fc £ fc (s) for s<Al (10) 

k=0 

where the coefficient is determined by the first k derivatives of U(s) at [i E . 
In our analysis the coefficient are all determined from ^ by the masses 
of the pseudo-scalar and vector mesons and the two parameters / and hp. 
It is emphasized that given our construction any contribution, be it a tree- 
level or loop term, contributes to the generalized potential exclusively via its 
derivatives at the expansion point. All left-hand cut structures are 'integrated' 
out systematically. In this work we construct the generalized potential U^ b (s) 
in terms of its first few derivatives at the expansion point fi 2 E , where we use 
the derivatives implied by the tree- level amplitude ([5]). Via (10) we obtain an 



approximate generalized potential for energies Ao < yfs < A s . For energies 
larger than the cutoff scale A s the generalized potentials are set to a constant 
[22] • By virtue of the specific form of the conformal map this is a smooth 
procedure. 

It remains to find the desired solution of ([7]). This is readily achieved in ap- 
plication of the N/D technique [30]. The partial- wave scattering amplitude is 
decomposed as 



T ab {s) = Y,D-J{s)N cb {s). (11) 

c 

Again, contributions of right- and left-hand cuts are separated. D ab (s) contains 
only the right-hand s-channel unitarity cuts. With the ansatz 



n /„n r s - /4f p(z>) v- r ds s-fi 2 M N ac (s)p cb (s) 



the coupled-channel unitarity is ensured [22]- In (12) we allow for one CDD 



pole structure [21] with its pole mass parameter M C dd and coupling matrix 
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. Such a term is requested by the presence of the s-channel vector meson 
exchange process in some J = 1 amplitudes |32,33j. The matrix N ab (s) con- 



tains left-hand cuts only. If the ansatz (11) is to represent a solution of the 



non-linear equation (|7| the matrix function N ab (s) has to satisfy the linear 
integral equation 
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The merit of the ansatz (13, 14) lies in its specification of the CDD pole 



parameters, R^ D ' and i?^, in terms of the parameters, g a ,Tn, characterizing 
a possible pole term in the generalized potential. They have to be dialed such 
that the effective potential U° b {s) is regular at s — m 2 . The CDD pole mass 
parameter, M C dd, is irrelevant. By construction, the scattering amplitude, 



which results from (11 14), does not depend on the choice of Mqdd- 



The system (13, 14) can be solved numerically by matrix inversion techniques. 



Once we obtain a solution of N ab (s), we can compute D ab (s) via (12), and a 
well defined result for the partial-wave scattering amplitude is obtained with 



(11). 



4 Numerical results 



Our results rely on the choice of various parameters. We use / = 90 MeV in 
the following, where we checked that small variations around that value lead to 
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Fig. 1. p-Wave scattering phases 5 for -kit — > irir with 1 = 1 and ttK — > irK with 
I = |. The data are taken form |34l35|36|37j . 

very similar results. For the cutoff parameter A s introduced in the conformal 
variable ^ we use the A s = 1.4 GeV unless stated otherwise. The small 
effects caused by variations around that central value will be discussed. The 
value of A s sets the scale from where on s-channel physics is integrated out. In 
our approach the generalized potential is characterized further by the number 



of terms consider in the conformal expansion (10). Typically a few terms 
only suffice to arrive at a good approximation of the generalized potential. A 
large number of terms would be justified if the left-hand cut structures are 
modeled very accurately, or a large number of counter terms are considered. 
If the empirical data set would be very accurate and more complete one could 



determine the various coefficients in (10) directly from the data set in a model 



independent manner. In our approach the left-hand cut structures are defined 
by the t- or u-channel exchange of Goldstone boson pairs. The latter imply 
specific correlations of the expansion coefficients. 

We begin with a presentation of partial- waves with J = 1, which probe the 
masses and decay properties of the light vector mesons. At tree-level the vector 
meson masses are 
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(15) 



which leads to the estimate mi- ~ 0.76 GeV and bo — 0.95. The latter values 
reproduce the empirical vector meson masses with an uncertainty of less than 
10 MeV. Within our unitarization scheme we identify the expressions of (15) 
with the mass parameter m in (14). The physical u and meson masses 
can be extracted from the (I G ,S) = (0~,0) sector, which involves the single 
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Fig. 2. s-Wave scattering phases for 
irK — > ttK with I = | (upper 
panel) and 7T7r — > tttt in I = 2 
(lower panel). The dotted, dash- 
dotted, dashed and solid lines cor- 
respond to a truncation at order 
0, 1, 2, 3 in the expansion (10) re- 



spectively. The data are taken from 
[36ll37l38l.'ffll4()l41l42l4.Slg] . 



channel KK only. The solution of (11 14) recovers the tree- level mass of the cu 



meson identically. Since its mass is below the KK threshold, its decay width 
is zero here. In contrast the (f) meson receives a small width depending on 
the actual value of hp. Its mass is almost indistinguishable from its tree-level 



expression (15). We determine hp ~ 0.26 as to obtain an accurate description 
of the (I,S) = (1,0) and (/, S) = (1/2,1) sectors, which are dominated by 
poles caused by the p and K* mesons. We find that the mass parameters m in 



( 14 ) needed to recover the empirical 7T7T and wK phase shifts are the empirical 
masses of the p and K* mesons. We conclude that the tree-level expressions for 
the vector meson masses survive the unitarization. This justifies our procedure 
to use physical vector meson masses in the driving expressions @. In Fig.l 
we confront our results with empirical p-wave tttt and ttK phase shifts. A 
good description up to about 1 GeV is achieved. The p and K* mesons are 
clearly visible in the data and our theoretical curves. As a consistency check 
we computed with ~ 3.47 MeV the width of the <ft meson. Our value is 
consistent with the empirical range of the partial width Y^KK) = 3.54 ±0.04 
MeV. It is amusing to observe that a universal value of hp is able to recover 
the decay properties of the p, K* and 4> mesons simultaneously. 

Our results are stable against variations of A s and changing the number of 



terms kept in (10). A variation of A s from 1.1 GeV to 1.6 GeV would almost 



not be visible in the figure and therefore we refrain from showing it. In contrast 
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Fig. 3. s-Wave scattering phases for 
tttt —J- 7T7T and tttt — > KK with 
1 = (upper panel) and the in- 
elasticity 7T7T — > KK (lower panel). 
The solid lines correspond to trun- 
cation the expansion 
The 



(10) 



terms. The data are 
jl5llH4l46l47l48l49l4.3l5()|51j 



with 4 
taken from 



s 1/2 [GeV] 



decreasing the number of terms in (10) from four to one defines the bands 
shown in Fig.l. 

We turn to a discussion of the s-waves, where we again focus on phases for 
which there is empirical information available. Our results for the 7r.fr phases 
are shown against empirical data in Fig.2. The figure includes also our result 
for the isospin two s-wave tttt scattering phase. While a variation of A s from 1.1 
GeV to 1.6 GeV would barely be visible in the figure, the results are sensitive 



to the numbers of terms considered in (10). The figure shows the results of 



keeping 1 to 4 terms in (10). In all three cases there is a saturation observed, 
with already the 4-term approximation recovering the empirical phases within 
the uncertainties of the different experimental results. We conclude from this 
result that indeed the vector meson exchanges contributions provide a quite 
accurate approximation for the leading coefficients in the conformal expansion 



(10). 



We close the presentation with a discussion of the isoscalar sector. Three differ- 
ent channels, tttt, K K and 7777 contribute in our analysis. In Fig.3 the empir- 
ical data from tttt scattering are collected and confronted against our results. 
The solid lines follow in the 4-term approximation. We recover the 7r tt — > 71 tt 
and 71 71 — > K K phases qualitatively, with some quantitative discrepancies 
around the / (980) resonance structure. The inelasticity in the tttt — > KK 
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reaction appears to be underestimated. 



Our results are in reasonable agreement with previous results from the Roy- 
Steiner equations [ofo]!?] . With our analysis we further explored the role of 
explicit vector-meson degrees of freedom. Our leading order computation has 
basically only two free parameters / = 90 MeV and hp = 0.26. To this 
extend we find the agreement with the previous results based on Roy-Steiner 
equations |5|6|7] the very satisfactory. 
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